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We derive forward inclusive dijet production in the scattering of a dilute hadron off an arbitrary 
dense target, whose partons with small fraction of momentum x are described by a Color Glass 
Condensate. Both multiple scattering and non-linear QCD evolution at small— a; are included. This 
is of relevance for measurements of two-particle correlations in the proton direction of proton-nucleus 
collisions at RHIC and LHC energies. The azimuthal angle distribution is peaked back to back and 
broadens as the momenta of the measured particles gets closer to the saturation scale. 



I. INTRODUCTION 



Understanding particle correlations in nucleus-nucleus collisions has been the purpose of many studies in recent 
\ years at the Relativistic Heavy Ion Collider (RHIC). These measurements provide insight on the properties of dense 
f^i . QCD matter, in particular they could help distinguish to what extend the experimental observations are due to 
initial-state or final-state effects. For instance, parton saturation [l[ is a characteristic of the initial nuclear wave 
functions, and is usually entangled with the final-state phenomenon of medium-induced energy loss Q. It is so in the 
suppression of high-p^ forward hadronic yields. By contrast, with measurements of particle correlations, initial-state 
and final-state effects may be disentangled. 

As a complement to single inclusive particle spectra, or to nuclear modification factors, the simplest observables are 
obtained with two-particle correlations. In particular, the correlation in azimuthal angle can provide further insight 
on the dynamics of energy loss. For instance it has been argued [1[ that the energy radiated off a hard final-state 
parton is expected to propagate in such a way that it should lead to an azimuthal correlation featuring a double-peak 
structure, instead of the standard back-to-back peak. This has been observed in Au-Au collisions, the away-side 
azimuthal correlation splits into a double peak with maxima displaced away from 180°. Meanwhile, the contribution 
of initial state saturation effects to the azimuthal decorrelation has yet to be determined. 

The main purpose of the present work is to study the inclusive two-particle spectrum for the following process: 
hT — > hih^X where h is a dilute hadron and the measured particles hi and /12 are detected in the very forward 
direction of that hadron. In this case, only the high-momentum valence partons of h contribute to the scattering 
and the dominant partonic subprocess is qT^qgX : a valence quark emits a virtual gluon g via lowest-order pQCD 
Bremsstrahlung and the quark-gluon fluctuation is put on shell by the interaction with the target T. In practice, this 
is relevant for forward particle production in deuteron-gold (d— Au) collisions at RHIC energies, where only valence 
, quark contribute. For proton-lead (p—Pb) collisions at the LHC, even in the forward region, it is likely that one also 
1 needs to account for the g^qq and <?— ><?<7 subprocesses. 

Meanwhile, it is suited to describe the target T by a Color Glass Condensate (CGC), as the process is mainly 
probing partons with a very small fraction of momentum x : the energetic probes formed by the projectile's valence 
quark, and quark-gluon fluctuation, interact coherently over the whole longitudinal extension of the target, and see 
the small— a; gluons inside T as a dense system of weakly interacting gluons, that behave collectivelly. It has recently 
become clear, and it is manifest in the process of interest here, that the CGC cannot be described in terms of a single 
. i-H , gluon distribution function, but rather must be described by n-point functions of Wilson-line operators. In other 
^ ■ words, the so-called fc-r— factorization used in [4] is not applicable. We explicitely derive the two-particle spectrum 
and obtain that it involves up to a 6-point function of Wilson lines in the fundamental representation. In addition, 
. . 1 we do not restraint the calculation to the soft-gluon approximation as done in [fj . 

We obtain a formula similar to that of [6(, but their counterpart of our n-point functions are treated differently 
than in the present paper. In this work, both multiple scattering and non-linear QCD evolution at small— a; are 
included: the Wilson- line correlators are computed in the framework of the Balitsky-Kovchegov (BK) equation 0]. 
Considering the correlations in azimuthal angle, we obtain that the perturbative back-to-back peak of the azimuthal 
angle distribution, which we recover for very large momenta of the measured particles, is reduced by initial state 
saturation effects. As the momenta decrease closer to the saturation scale Q s , the angular distribution broadens, but 
at RHIC and even LHC energies, saturation does not lead to a complete disappearance of the back-to-back peak. 
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The plan of the paper is as follows. In Section II, the inclusive two-particle spectrum qT^qgX is calculated. The 
cross-section, differential with respect to the quark and gluon transverse momenta and rapidities, is expressed in terms 
of Wilson-line correlators in the CGC wavefunction. Section III is devoted to those CGC averages and describes how 
to perform them in the context of the BK evolution. In Section IV, we introduce the relevant observable to study the 
azimuthal angle correlation and we investigate the impact of the small— a; evolution at RHIC and LHC energies. 



II. FORWARD INCLUSIVE DUET PRODUCTION CROSS-SECTION 



In this section we derive the inclusive quark-gluon production cross-section in the high-energy scattering of a quark 
off a Color Glass Condensate. We shall work at leading order with respect to the strong coupling constant as, and 
to all orders with respect to gsA where A is the CGC color field. The relevant diagrams are shown in Fig.l and the 
necessity of both diagrams will become transparent in the following derivation: it is a manifestation that there should 
be no gluon production without interaction with the target. We will later require that the scattering process features 
a hard momentum transfer A> A-qcd, which justifies the use of perturbation theory. 

We shall use light-cone coordinates with the incoming quark moving along the x + direction, and work in the light- 
cone gauge A + = 0. In this case, when the quark passes through the CGC and interacts with its color field, the 
dominant couplings are eikonal: the partonic components of the dressed quark wavefunction have frozen transverse 
coordinates and the gluonic field of the target does not vary during the interaction. This is justified since the incident 
projectile propagates at nearly the speed of light and its time of propagation through the target is shorter than the 
natural time scale on which the target fields vary. The effect of the interaction with the target is that the components 
of the dressed quark wavefunction pick up eikonal phases. 



A. The dressed quark wavefunction 

To describe the wavefunction of the incoming quark, we shall use light-cone perturbation theory. We denote the 
3— momentum of the quark p = (p + ,p±), and its spin and color indices a and i. One has p~ = (p\ + m 2 )/2p + . To 
decompose its wavefunction on bare quark and gluon states, we introduce b ia (k) and bi. a (k) (resp. a^ A (fc) and 
Qc,\(k)), the creation and annihilation operators of a quark with color i, spin a (resp. gluon with color c, polarization 
A), and 3— momentum k. One has 

b\jk)\0) = \k,i,a) , Mfc)|0) = , {M*0, &],,#')} = <M^ (3) ( fc ~ k ') , W 
4,a(*)|0) = \k,c,X) , a C)A (fc)|0) = , [a c ,x(k), 4,v( fc ')] = S cd 6 xy 5^(k - k') . (2) 

We recall that in light-cone quantization, the virtuality of the particle comes from the non-convervation of the 
momentum in the x~ direction, meaning (k— k')~ =/=k~ — k'~ . All particles are on-shell and only the 3— momentum is 
conserved. 

We work at lowest order with respect to as so we only need to consider the Fock states \q)o and \qg)o hi the 
decomposition of the dressed quark: 

\p,i,a) = Z\p,i,a) Q + d 3 k g s T^ ip^pip, k) \{p-k, j, (3); (k, c, A)) . (3) 

The \qg)o part of the dressed quark is characterized by the wavefunction gsTfjipap where (3 and j denote the spin and 
color indices of the quark and A and c denote the polarization and color indices of the gluon. T c is the generator of 
the fundamental representation of SU(N C ) and ip^p 1S gi ven by 

1 upip-k^^A^u^p) 
^ k) = V8(p-fc)W {p-k)- + k-- P - (4) 

where k = (k + , k±) and p—k = (p + — k + ,p± — kj_) are the 3— momenta of the gluon and quark respectively. The factor 
Z is a renormalization constant determined from the requirement that the normalization of the dressed quark is the 
same than that of the bare quark (note that those normalizations are proportional to 5^ (0) because we are using 
plane waves). Physically, Z accounts for the virtual gluon emission associated with the real gluon emission explicit 
in When calculating inelastic processes, as is the case here, the actual value of Z is not relevant. 



FIG. 1: Inclusive quark-gluon production cross-section in the high-energy scattering of a quark off a Color Glass Condensate. 
p = (p + ,p±) : momentum of the incoming quark; q= (q + , q±) and k = (k + , k±) : momentum of the outgoing quark and gluon. 
The vertical wavy lines represent the interaction with the target CGC, each line carries a factor gsA and multiple gluon 
exchanges must be resummed. The black points represent the emission of the produced gluon by the quark, it is emitted before 
the interaction or after the interaction in which case the contribution comes with a minus sign, as explained in the text. 



With our choice of gauge, one can write the two gluon polarization vectors e^. (k) and e^s (fc) in terms of two 
transverse vectors e]_ and e\ : s^Jk) = (0,6*1, k± ' £ ±A' + )- Then using the chiral representation of the Dirac 
matrices to obtains the spinors leads to 



V2(kj_-zp±) ■ e 1 L [5 a -5 /3 - + (l-z)5 a+ 6p + ]+m,z 2 6 a+ Sp- A = 1 
V2(kj_ — zp±) ■ e 2 L [5 a+ 5 f 3 + + (l-z)S a -5 l 3-]-mz 2 6 a -5fj + A = 2 



(5) 



with z = k + /p + . For reasons which will become clear later on, it is more convenient to work in a mixed space, in 
which the transverse momenta are Fourier transformed into transverse coordinates: 



b{ a (p+,h) = J d 2 b e-^ b &J ia (p+ j>j.) , < A (p+ x) = / d 2 b e-*^ 



PA 



(6) 



In the mixed representation, the decomposition of the dressed quark ^ on the Fock states |g)o and \qg)o is the 
following: 



\p,i,a) 



d 2 b 



Z\p+,b,i,a) +J2 I * + 7TT2 9sT? ^(p,fc + ,x-b) \(p+-k+,b, 



3, 13); (fc + ,x,c, A)) 



(7) 

where b and x are the transverse positions of the quark and gluon respectively, and with the mixed-space wavefunction 



^ /3 (p,fc+,x)= / d 2 k ± e lk ^^(p,k) 



(8) 



given by 



X • £ | 



b* (p,k+,x) = 



<k+ 



|x| 

izV2Ki(mz\x.\)— —^[5 a +5{3+ + (l — z)5 a -8p-] — z 2 K (mz\x.\)S a -Sf3+ A = 2 

3C 



{6 a -6f3- + (l-z)6 a+ 6p + }+z 2 K (mz\x.\)6 a+ 6i3- A = 1 



(9) 



B. High-energy eikonal scattering off the target 



Let us first recall the basics of the CGC description. When probing inside a target hadron with processes that 
are sentitive to partons with a very small fraction of momentum x, the probe actually sees a dense system of gluons, 
responsible for large classical color fields A ~ 1 / gs [8] . Rather than using a Fock-state decomposition which is not 
adapted to account for the collective behavior of the small— a: gluons, it is more appropriate to use other degrees of 
freedom and describe the target by classical color fields: 



\T) = \qqq)o + \qqqg)o + ■■•+ \qqqg ■ ■ ■ ggg)o + ■■■ => \T) = DA $ XA [A] \A) . (10) 



The CGC wavefunction $^^[^4] is normalized such that J DA\<f> XA [A]\ = 1, and xa denotes the smallest fraction 
of longitudinal momentum probed. It depends on the final-state kinematics of the process condidered and will be 
specified later. The CGC wavefunction is mainly a non-perturbative quantity, but the x a evolution of |4>a; A [*4]| 2 can 
be computed from perturbative QCD [9(, in the leading— 1u(1/xa) approximation that resums powers of as Iii(1/xa)- 
It is a priori not obvious that this description of the target, which requires small values of xa, is valid for experiments 
at present energies. However, it has had success for many observables in the context of HERA [lGj and RHIC [111 ]. 

The target is moving along the light-cone in the x~ direction, and the only component of its color current is J~. 
With our choice of gauge A + = 0, the current conservation law [D^, J M ] =<9 + J~ =0 implies that J~ does not depend 
oni". Therefore one writes 

J"(x") = 5»-J-(x+,x) = 5»-T c p c (x+,x) (11) 

where we have introduced the color charge density of the target p c . Solving the Yang- Mills equations [Dp, F^} = J v 
then leads to (see for instance 

A*(x u ) = 8^-T c A-{x+,x) , -V 2 A-(x+,x) = p c (x+,x) . (12) 

The formal functional integration in (|10[) stands for the A~ integration. 

In a scattering process, the outgoing state is obtained from the incoming state by action of the scattering matrix 
S. When high-energy partons scatter off the CGC, the interaction is eikonal and the S matrix acts on quarks and 
gluons as (see for example [H[): 

S\b,i)®\A) = '$2w$[A]Qa)\b,j)®\A) , S\x,c)®\A) =Y,W A d lA}(x)\x,d)(g>\A) , (13) 

j d 

where the phase shifts due to the interaction are described by Wf and Wa, the eikonal Wilson lines in the fundamental 
and adjoint representations respectively, corresponding to propagating quarks and gluons. They are given by 



W f ,a[A](x) = Vexp (ig s J dx+T F A A^ (x + ,x) 



(14) 



with Tp A the generators of SU (N c ) in the fundamental (F) or adjoint (A) representations and with V denoting an 
ordering in x + . The Wilson lines resum powers of gsA which is necessary as A is a large classical field whose strength 
is of order l/gs- It is manifest from (| 13|> why the mixed representation introduced earlier is convenient: working with 
transverse space coordinates (instead of transverse momenta) provides eigenstates of the high-energy S— matrix. 
Coming back to our computation, the incoming and outgoing states are 

\V in ) = \p, i, a) ® \T) , \* ut) = 5|* in ) . (15) 

Using (|7|) and p3p . it is straightforward to obtain (we now keep the A dependence of the Wilson lines Wf.a implicit) 

!*<>«*> =JdA $ xa [A J ^ e^' b (z^\W F (b)] ij \p + ,bJ,a) Q ®\A)+ 
V f dk+^ as6*Jp,k + ,x-b)[T c WF(h)hW? 



gs^(p,k + ,x~b)[T c W F (b)] t3 W^ d (x) |(p + -fc + ,b,j,/3);(fc+,x,d,A))o® \A) . (16) 



In this formula, the quark-gluon contribution represents the first contribution pictured in Fig.l, for which the gluon 
is emitted before the interaction. The second contribution, for which the gluon is emitted after the interaction, is 
hidden in the quark contribution Z\p + , b, i, a)o- To see that, let us write this in terms of \p,i,a) : from formula [|7]). 
and using the fact that the p± dependence of <^g(p, k + , x) is exp (ik + p±.x/p + ), one obtains: 

Z\p+,b,i,a} = f d 2 p ± e-^- h \p,i,a) - f dk+^d 2 b' d(b-b' ^k+(x~b')/p+)g s T^ e ^-Cb'-b) 

^(p.A+.x-b') |(p + -fc+,b',j,/?);(fc+,x,c,A)) • (17) 

One sees that the emission-after- interaction term arises with a minus sign. The dressed quark contribution \p, i, a) 
does not contribute to gluon production, and can be removed from the final answer. Indeed, when computing gluon 
production, gluons which dress the final-state quark should not be included, as these are not actually produced. 



Combining (|16p and (117|) , the outgoing state can be simply rewritten as 

|*«*>= j DA® XA [A\Y. J dk+ ^]0j2 eiP± ' h 9s K X ^ k+ ^ h ) l(p + -fc + ,b,j,/3);(fc+,x,c,A))®|^) (18) 
with 

^ )i >,fe+,x,b) = ^(p.fc+.x-b) [T d W/ F (b)t<(x) - W/ F (b + fc+(x-b)/p+)T c ] y . . (19) 

Once again, the different contributions contained in this wavefunction have a straightforward physical meaning: the 
term containing W a corresponds to the interation happening after the gluon emission while the contribution without 
Wa corresponds to the interation taking place before. 



qT ^qgX 



C. The inclusive quark-gluon production cross-section a 

From the outgoing state (|18|) . one can now compute the production of a quark with momentum q = (q + , q±) and a 
gluon with momentum k — (k + , k±). The corresponding cross-section cross-section reads 



Lj2(* out \N q (q)N g (k)\y out ) (20) 



da qT^qgX 

d 3 kd 3 q ~ 2N t 

where we recall that i and a refer to the color and spin of the incoming quark of momentum p = (p + ,p±). In (|20l) . 
the operators N q and N g are given by 

N q {q) = E b lA^Al) , N g (k) = Y, ol, x (k)a a>x (k) (21) 

in terms of the creation and annihilation operators introduced earlier in Section II-A. Let us rewrite the cross-section 
(|20p using operators in the mixed representation 

'/"" J -"' A i r >i 2 ->-' <?b d%> eifcx . (X '_ x)eigx . (b '_ b) 



d 3 kd 3 q 2N C J (2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 

E (* out |ai,A(x',fc + )^ /3 ( b '^ + ) a -A(x,fc + )^, /3 (b,g + )|*out) , (22) 

aAa/3zj 

where x and x' (resp. b and b') represent now the transverse coordinates of the measured gluon (resp. quark) in the 
amplitude and the complex conjugate amplitude respectively. When computing the action of a a ^\ and bj.p on \^ ou t), 
one obtains 

a ,A(x,fe+)6 j!/3 (b,g+)|* 0Ut ) = g s e^- b S(p+-k+-q+) J DA $ XA [A] ^,y(p, k+, x, b)\A) . (23) 

The function 6(p + — k + — q + ) in (|2"3"|) shows that the longitudinal momenta are converved during in the high-energy 
eikonal scattering. However, this delta function at the level of the amplitude will lead to a factor 6(0) when computing 
the cross-section (f2"2"j). This problem is related to the factor 6( 3 )(0) present in the normalization of the state l^m), 
itself due to the fact that we are using plane waves. When computing a physical observable, this divergence usually 
goes away, as is the case for the two transverse dimensions. However, a 6(0) remains for the longitudinal direction 
precisely because longitudinal momenta are conserved by the interaction. Working with wave packets would solve the 
problem, and the appropriate prescription is to remove the factor 2ir6(p + — k + — q + ) from the amplitude (|23p . and to 
put it back in the cross-section (f2"2"|) . 

Denoting z = k + /p + and using (f2"3"f in (f2"2"|) . one obtains the qT^qgX cross-section (as = g|/47r): 

d ^ 9X -a s C F 6(p+-k+-q+) e ^, X '- xV( ^), b '- b ) 



d 3 kd 3 q -a"*"**' * J J ( 27r )2 (2tt) 2 (2tt) 2 (2tt) 

E^fe k+ -h')^^, fc+,x-b) |s*^ g [b,x, b',x.';x A ] - 5^[b,x, b' + z(x'-b'); x A ] 



Aq/3 



■S2[b+ Z (x-b),x',b';i A ] + s(f ) [b+z(x-b),b'+, 2 (x'-b');z A ]} • (24) 



In we have introduced the following traces of products of Wilson lines: 



^)_ ff ( b ,x,b',x';zA) - — (Tr(iy F (b)^(b')T d T c )[T^ A (x)^(x')] c<i )^ , (25) 
S^(b,x,b';^) = ^-L- (Tr(^(b')T c ^(b)T d )^(x))^ , (26) 

S%(b,b';x A ) = i- (Tr(^(b)^(b')))^ , (27) 
and we have also denoted the average over the CGC wavefunction squared |$a; A [./4]| 2 using the following notation: 

J DA\$ XA [A]\ 2 f[A] = {f) XA . (28) 

The quantities Sj^ , and 5*^^ g co ntain the QCD evolution toward small values of Xa- 
Some comments about formula (124p are in order. 

• First one recovers the result of except that in our case the interaction with the target is obtained by 
specific n— point functions, expressed in terms of Wilson lines that resum powers of gs-A. These are to be 
averaged with the CGC wavefunction squared [„4]| 2 , whose xa evolution resums powers of as ln(l/x A ), in 
the leading— ln(l/a;^) approximation. 



• Using the Fierz identities 



[WHx)] y [Wl( x )] kl = i- 6 a S jk + 2W c A d (x)T^ , (29) 

T ij T u = \^i 5 3k - yySijdki , (30) 

one obtains that W A d (x) = 2Tr(JU F (x)T c WF(x)T d ) which shows that an adjoint Wilson line Wa is equivalent to 
two fundamental Wilson lines Wp. Therefore the quantities sjjg , ^q 3 g\ ancl Sqg\g are ^~ ana - 6-point functions 
with respect to the averaging l|28p. Using (f2"9"|) and ([3"U|) . one can also see that the singularities of (f>*g in ([2"4"]) 
are cancelled when x = b or x' = b' . 

• Our result should be identical to that of flij . although this is not so straightforward to see, as their expression 
for the qT — > qgX cross-section is not as compact as our formula ([24)l . This is probably because in [lj], the 
Wilson lines are expressed in momentum space, where the interaction is not diagonal. 

Let us finally consider the soft-gluon approximation 2<1. Not only does the wavefunction ^ simplify, but also 
many simplifications occur with the Wilson lines. Indeed using the Fierz identity to transform all the adjoint Wilson 
lines into fundamental ones, one obtains (if one also factorizes the remaining traces, meaning (Tr( . )Tr( . )) XA = 
(Tr( . )) XA (Ti( . )) XA , we recover the expression of [5j): 



Sl%(b,x,b'y;x A ) - sl%(b,x,b';x A ) - S q %(b, x', b'; x A ) + sg>(b, b'; x A ) = ^ /^Tr(V F (b)IU F (b') 



J-Tr(w F (b)Wt(x))i-Tr(^(x)^(b')) - ±-Tr(w F (b)W F (x')) i-Tr(V F (x')W F (b')) 

+ ±Tr(lU F (x)IU F (x'))^Tr(lU F (b)IU F (b')IU F (x')IU f (x))\ .(31) 



N, 



The last term in this expression shows that even in the soft-gluon approximation, the dipole degrees of freedom 
(i.e. traces of two Wilson lines) are not sufficient to compute the dijet cross-section ((24"|) . This is an important 
difference with the case of single-particle production for which using dipoles is enough, yielding the so-called 
kx— factorization. For two-particle production, kx— factorization cannot be used. 

In the following we shall not consider the soft-gluon approximation, and work directly with the 2-, 4- and 6-point 
functions ([2T|) . (f2"6"|) and ([^5]) . Indeed, we are interested in final-state configurations where the two particles are 
both detected in the forward hemisphere, with similar rapiditites. By contrast, using a soft gluon would impose the 
constraint of a large rapidity between the two detected hadrons. 



III. PERFORMING THE TARGET AVERAGES 



In this section, we explain how to model the CGC wavefunction in order to perform the average while taking 

into account the small— a; QCD evolution. To compute the correlators and S gg \ g , we shall model the CGC 

wavefunction squared [*4]| 2 using Gaussian-distributed sources. Then, following the approach of [lq . we will work 
in the large— -/V c limit which allows to easily implement the Balitsky-Kovchegov (BK) evolution, and also simplifies 
the analytic expression for the 6-point function. 



A. A Gaussian distribution of sources 



To compute the average (|28|) . we use the following Gaussian distribution of sources 



where the color charge density p c and the color field A are simply related via formula p2|) . The variance \i XA is a 
function of xa and caracterizes the density of the color charges. With this Gaussian approximation, the 2-point and 
4-point correlators (|27|) and (|26[) have been computed for arbitrary N c (see for instance the Appendix A of They 
are given by 

sif(b,b';x A ) = e-^rcb-b',,.) ) (33) 
sg(b,x,b';xA) = e -^(x-^A)+r(x~ b ', XA)]+jk r( h ~ h ', XA) ^ (34) 

with the function r(b — h ,xa) related to fi^ A in the following way 

r(b-b',au)=4 [d 2 xd%dz+ »l A (z + , X -y)[G (b-x)-G (b'-x)][G (y-b)-G (y-b')} , (35) 



where Go is the two-dimensional massless propagator 

d 2 k e lkx 
(2^T 



Go(x)= / j^-^. (36) 



|k|>Aq C D 

The last correlator needed for our study has a more complicated structure. Following the derivation of the 4-point 
function in [13], it is easy to see that the problem of computating the six point function ([25| can be reduced to the 
diagonalisation of a 6 x 6 matrix. Note that in the recent publication jl8| . a method to compute any n— point function 
is presented. As already mentioned, only the large— N c limit is used in the implementation of the small— X QCD 
evolution. Therefore, computing the exact structure of the 6-point function S^q g is of no interest for the present 
study. We only give the large— N c result 

Sg ) ffl (b,x,b',x'; a u) =e -^[r(x-b,x A )+r ( x'-b',x A )+r(x- x ',^)] _ (37) 

This will be used along with formula ((3"3"|) and the large— N c limit of formula (j3~4")l . Note that these expression should 
be understood for scattering at fixed impact parameter, which is why they feature one less independent variable than 
excepted. We will assume in the following that the impact-parameter dependence of the correlators factorizes. 



B. Evolving the MV model with the BK equation 

Let us now explain the strategy to evaluate the function T(r, xa), in terms of which all the correlators (|3"3"|) . (f3"4")) and 
(l3"Tf could be written. We will use the fact that, under the asumption |3"^|) that the CGC wavefunction is Gaussian, 
the large— N c limit implies the following result 

(lt(w^(b)wJ(b / ))lt(w>(x)W+(x / ))) =( I &(w F (b)T4(b'))) (t>(wHx)W£(30)) • (38) 



This significantly simplifies the high-energy QCD evolution equations. Indeed, when considering the Balitsky hierarchy 
of equations for the n— point correlators, (which is a rewritting the JIMWLK functional equation for the evolution 
of \&x A [-4] | 2 with xa), the factorization ([38]) reduces the QCD evolution to a single closed non-linear equation for 



Sgg (b, b'; Xa), known as the BK equation [7|- In this work, we consider the impact-parameter independent version 



j(2) 

dS™ (b-b';x) _ f tPz (b-b') 2 / (2) z . x)s (2) ( h , ) s (2) (h b , A m) 
d\n(l/x) ~ a J 2tt (b-z)2( z -b') 2 V^ (b Z '^^ Z b ' X > ^ b b ' X >) ^ 

with a = asN c /ir. From formula (135|) . the solution of the BK equation (|3^|) gives L(b — b ,Xa), which allows to 
compute all the n— point functions needed for the calculation of the qT^qgX cross-section. 

As the initial condition, we shall use the McLerran-Venugopalan (MV) model Q for T(r, xq). It is obtained from a 
Gaussian average of the type (j3"2")l , with /j, 2 , o (z + ,x) = £(x)/i 2 o (z + ) : 

T(r,x ) = fir* (J dz+ M 2 o (z+)) £ dk . (40) 

For |r|AQCT> "C 1, the remaining integral behaves as — ln(r 2 Ag C£) )/(87r). This leads to the MV model: 

r Mo) = _i_ r * e ^ e + _l_), Ql = ^j dz+<(z ^ (41) 

where Q So is the initial saturation scale, at x = xq. In practice, we choose to start the BK evolution at xq = 0.01, which 
is small enough to justify using the MV model. For a target Gold or Lead nucleus, 2-kQ 2 = 2 GeV 2 is appropriate 
at xo = 0.01 (we use the same value as in [l6j], our saturation scales differ by a factor V2tt). For that moderately 
small value of x, the MV model gives reasonable results in heavy ion collisions at RHIC [HI]. For x<xq, the quantum 
evolution effects are then implemented by the BK evolution. 

It is known that the BK equation (|3U)) . equiped with the leading-logarithmic Balitsky- Fadin-Kuraev-Lipatov 
(BFKL) kernel x(l) = 2^(1)- VK7) — V^l - 7) (i n Mellin space), leads to an increase of the saturation scale that goes as 

Q» = Ql (-) Va with v = x'hc) = ^ = 4.88 ( 7c = 0.6275) . (42) 



' S0 \x; lc 



In practice however, the appropriate saturation exponent is Q 2 s ^x~ x with A~0.25, and this discrepency is understood 
in terms of subleading logarithms [20] | . In this work, in order to mimic to correct evolution of the saturation scale 
with x, we impose the unphysical value a = 0.05 when solving the BK equation (|39[) . This is a simple way to account 
for the unknown next-leading effects, while staying compatible with experimental observations. 



C. Expression for the cross-section 

We now come back to the expression (|24|) for the qT^qgX cross-section, which can be simplified using the results 
([531) . (GS) and dSZJ) for the 2-, 4- and 6-point functions (gTJ), and ([23 in the large-iV c limit. Let us first factor 
out some prefactors: 

da qT ^ qgX 

d3M3q =a s C F S(p + -k+-q + )M(p,k,q) . (43) 
Changing the integration variables to u = x b, v = zx+(l — z)h, u'=x' — b , and v' = zx' + (l — z)h , one writes 

"wi-/w' t ' H '5* k,t ^ ,t ' ,1 /ffi i P^ 

^ e -^[T(u,a; J i)+r(u',a; A )+r(v-V+(l-a)(u-u'),aiA)] _ e -^[r(u,x A )+r(v-v'+(l-2)u,^)] 

_ e -^[r(.M',x A )+r( v '- v +(i-z)u\x A )] + e -i^r( v - v \x A )^ ^ 

where we have introduced the following transerve momenta 

K = (1 - z )k± + z(p± - q±) , A = k± + q± - p± . (45) 



In (|44|) , the transverse momentum transferee! during the collision A is Fourier conjugate to v— v', while k, which 
caracterizes the invariant mass of the final-state qg system, is Fourier conjugate to u— u'. 

We now change the integration variables v and v' to r = v — v'. and B = (v + v')/2. The integration over B 
represents the impact parameter integration. Following our approximation that in the correlators the B dependence 
factorizes (and is not explicitcly indicated), one has J d 2 B(l — S) = St (1 — 5) and the B integration simply yields 
the normalization factor St, which caracterizes the transverse area of the target. One obtains 



Ofpi -5i)-(u'-u) g - ^ [r(u,a;A)+r(u',i j)] _ e iK-u' -i(p ± -q ± )-u e -^[r(u,x A )] 

-iK-u+i(j>j_-qj_)-u' — ^[r(u',£CA)] _l „i«-(u'-u 



_ e -lK,-u+t(p ± -q±)-u g-^-l 1 \ u > X A)\ _|_ g«-(,u -u) j _ ^g-j 

Finally, performing the Fourier transforms back to momentum space, the cross-section can be written in the following 
compact form: 

^ = ST ^ 6(p+ - k+ - q+) S \^ k + ,Pt-q^x A ) - V^(P, *)\ 2 F XA (A) (47) 

where we recall that is the q^qg wavefunction ([5|) in momentum space, and with 

I^(p, fc + , «; a*) = | 70^- ■ u e-^ r ("^ fc + , u) . (48) 

We have also introduced the so-called unintegrated gluon distribution F XA (A), which is simply the Fourier transform 

(2) 

of the 2-point function S q q : 

F XA (A) = l^e-^e-^^. (49) 

It is important to stress that F XA (A) is not the usual unintegrated distribution defined in the leading-twist ap- 
proximation of perturbative QCD, but it is rather an all- twist gluon distribution, as multiple scattering are taken 
into account. Although they coincide for large momentum, ^^(A) contains more information, as it is also properly 
defined in the infrared (meaning for A 2 < Q 2 S ). Note that the numerical Fourier transformation in (|4"9"f suffers from a 
positivity problem and features oscillations for large momenta [2lj | . This is due to the infrared behavior of the initial 
condition T(r, xq), and the problem is carried to lower values of xa by the BK evolution. To cure this, we substituted 
the 0(fc — |r|AQC'D) function in (|4"0")) by the smooth cutoff function 2 arctan[(/c/|r|AQc:z?) 2 ]/7i'- This modification of 
the infrared regularisation does not influence the results provided F XA is evaluated for values of |A| much bigger than 
Aqcdj and this will be the case 

Finally let us insist that, as already mentioned, the unintegrated gluon distribution F XA is not enough to characterize 
the CGC, and for the process considered here, more information is needed to compute the cross-section. It is contained 
in the modified wavefunction l\ : 

1 ( V2e]_ ■ G XA (K-zpj_)[6 a -6p- + (l-z)6 a+ 5p + ]+H XA (K- zp±)5 a+ 6p- A = l 
I^{p,k + > K;x A ) = -^i (50) 
V [ \[2e\ ■ G xa (k- 2p±)[6 a+ 6p + + (l — z)6 a -6p-] — H XA (K—zp±)6 a -5p+ A = 2 

given in terms of the functions 

G xa {k) = mz J^f du e-^ LT{u ' XA) uK 1 (mzu)J 1 (\K\u) , (51) 
H xa (k) = mz 2 { du e~^ T{u ' XA) uK a {mzu)Ju{\n\u) , (52) 



that contain the rest of the xa dependence. The functions F XA , G XA (which is a two-dimensional vector) and H XA 
are all obtained from the function Y(u= |u|, xa), itself obtained by solving the BK equation. 



IV. AZIMUTHAL ANGLE CORRELATIONS IN d-Au COLLISIONS 



In this final section we consider the inclusive two-particle spectrum for the process hT — > h\hiX and, as an 
application, we study the correlations in azimuthal angle between the measured particles h\ and /12, which are 
both detected at forward rapidities. We consider deuteron-gold collisions at RHIC center-of-mass energies (yfs = 
200 GeV/nucleon) and investigate the suppresion of the back-to-back peak as a function of the transverse momenta 
and rapidities of the two particles. 



A. The inclusive quark-gluon production cross-section o hT ~^ q9X 

Let us denote P + the momentum of the hadron h. To obtain the cross-section a hr ^ qgX from the partonic cross- 
section a qT ^ qgX , one can use the collinear factorization of the quark density inside the hadron q(p + / P + , /1 2 ) where 
/i 2 is the factorization scale: 

rf n hT^qgX f J qT^qgX 

dx q(x,fi 2 ) „,„ (p+=xP+,p ± =0) . (53) 



d 3 kd 3 q J ' d 3 kd 3 q 

Indeed, we are interested in measurements that probe only large values of p + /P + . By contrast, partons with small 
fraction of momentum are probed inside the target, and it cannot be described by a single gluon distribution (as 
obvious from formula (f2~j| for instance) . We shall later use the factorization scale /i 2 = A 2 which is the transverse 
momentum transfcrcd during the collision, and it is supposed to be large to justify our perturbative calculations. 
As a function of k± and A, the a hT ^ qgX cross-section reads 

= S T ^^(l-z)x h q(x h ,A 2 )k + T\l x f3 {p 1 k + ,k ± ~A;x A )~^ f} (p,k + ,k ± ~zA)\ F XA (A) (54) 



d 2 k ± d^dy k dy q -*r-Z£-V-*,Wi*~ J- 

where in I X p and P should be understood as (p + =q + +k + ,p± =0). Also, z now stands for z = k + /(k + +q + ) and 
we have denoted Xh = (k + +q + )/ P + , the fraction of momentum of the probed quark inside the incoming hadron. Note 
that similarly, there is conservation of momentum along the x~ direction and one has xa = {k~ +q~)/ P~ where P~ 
is the momentum of the incoming target. In terms of the rapidity of the quark y q , the rapidity of the gluon y q and 
the center of mass energy y/s, one has: 

\kx\e y " |fc±|e»* + |g±|e»« Ifcxle - ** + |gj_|e-»« 

Xh = 7= , x A = 7= ■ (55) 



\kA_\eV" + \q±\ev* ' 

It is clear that in order to have Xh < 1 and xa *C 1, one needs forward rapidities, in the hemisphere in which the 
hadron h fragments, where y q and yk are positive. 

In this work, we do not take into account the fragmentation of the final-state quark and gluon into hadrons. Since 
we are interested in azimuthal angle correlations, fragmentation does not play a important role. Also, we do not 
include in the process hT ^h\h%K, the contributions of the gluon-initiated subprocesses gT — > qqX and gT — > ggX. 
Measurements at forward rapidities at RHIC energies involve values of Xh so high that they are only sensitive to the 
valence quarks. In the case of the LHC, it is likely that one also needs to account for the gluon-initiated subprocesses. 

In order to simplify the numerical computations, we shall work with massless quarks. Using m — yields 

G XA (k±) = ^-7 [du e-^^JiClAjJu) , H XA (k±) = (56) 



\k 

and in (|54j) . the summation over the quark and gluon spins and polarizations becomes 



k + ^2\l^(p,k + ,k x -Aix A ) - ii X p(p,k+,k ± -zA)\ = 2 [l + (l-z) 2 ] 

Aq/3 



G XA (k ± -A) 



kt-zA :J 



|fc x -^A| 2 



(57) 



In Fig. 2, we display the functions F XA and G XA (in the massless case) obtained from the BK evolution. More precisely, 
the function A 2 F XA (A) is displayed in Fig. 2a and as is well-known, it is peaked for A 2 ~ Q 2 , and therefore as xa 
decreases the peak travels towards higher momenta. The evolution is however quite slow, considering the small value 
of the saturation exponent A. The function 1 — k± ■ G XA (k±) is represented in Fig. 2b, it decreases from 1 to as k\ 
increases, and the front travels to higher momenta as xa decreases. 
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FIG. 2: Left plot: the dimensionless unintegrated gluon distribution A 2 F XA (A) as a function of A 2 (see formula (|49p ). Right 
plot: the dimensionless quantity 1 — k± ■ G XA (k±) as a function of k\ (see formula (|56p . in the massless case). The curve for 
xa~ xo = 0.01 is obtained from the MV initial condition (|40p . The evolution towards smaller values of xa is obtained with the 
BK equation (|39[) . and is shown over 6 units of ln(a;o/a;A). 



At this point, it is easy to recover the perturbative limit of the inclusive two-particle spectrum for \k±\^>Q s , 
Gx A (k±_) = fc±/|fc±| 2 leads to k T — factorization. For |A|^>Q S one has F XA (A) = Q 2 (x a) / A A ) [22| which yields 

Klch 9X a s C F {l-zni + {l-zf]Ql{x A ) 

Ot—^t-s T, w> \?n XT? XhQ(Xh,A ) . (58) 



d 2 k ± d 2 q ± dy k dy q 2tt 3 {k A _-A) 2 A 2 (k ± -zAy 



B. An application: azimuthal angle decorrelations 

We will now use the inclusive two-particle spectrum to investigate the process hT^h\h 2 X, and in particular 
the cross-section as a function of A(f> = 4>\ — <\>ii the difference in azimuthal angles of the two measured particles h\ 
and h 2 . We will study the normalized A</> distribution 

1 da f da hT ^ h ^ x y 1 da hT ^ h ^ x 
a dAip ~ \ dp Tl dpr 2 dyidy 2 J dp Tl dp T2 dy x dy 2 dA<j> 

where {pT t , 4>i) an d {pt 2 > ^2) are the transverse momenta of the measured hadrons and y\ and y 2 are their rapidities. 
Our results can be applied to d—Au collisions at RHIC, and to compute Xh and xa we will use -^=200 GeV. 

A given final-state configuration can be obtained in two possible ways from the hT —> qgX process, depending on 
which particle (1 or 2) comes from the quark and which comes from the gluon. While Xh and x A are the same in both 
situations, it is not the case of z (which is changed into 1 — z), and therefore the cross-section (|54|) is not symmetric 
with respect to the two situations, as it is a decreasing function of z. Let us choose to label the particles such that 
PT^eV 1 >PT 2 eV2 - If the quark (resp. gluon) is the particle 1, then z< 1/2 (resp. z> 1/2). This shows that, for similar 
transverse momenta pt x r * J PT 2 i the favored configuration is the one where the quark is the most forward particle; this 
is especially true when y\ — y 2 is large. In any case, we take into account both situations. 

The massless quark approximation is valid when (fci-zA) 2 >m 2 , therefore we will stay away from the situation 
2/1=2/2, because when A(\> = it leads (k± — zA) 2 = 0. In this situation the factor (k± — zA)~ 2 in (|57|) should actually 
be replaced by 1/m 2 ; there is an increase of the cross-section when A0 ~ and y\ —y 2 , which corresponds to the 
quark and the gluon being collinear. Also we shall not consider the situation px l =Pt 2 which implies that A 2 = for 
A(/) = tt. Indeed, we would like to work with |A|^>Aqct>- 

As can be seen from the kinematics the most forward of the two particles essentially determines the value 

of Xh while the most central one determines the value of xa- In order to study the effect of the CGC evolution, the 
ideal situation would be to keep Xh fixed and to vary xa ■ In practice, this is better realized by fixing the rapidity and 
momentum of the most forward particle and by varying the kinematics of the other. Note that doing the opposite 
would emphasize the Xh evolution of q(xh,A 2 ), rather than focus on the xa evolution of F XA and G XA . Moreover, 
the cross-section ()54[) is quite sensitive to choice of factorization scale in the quark density, so it is better to keep Xh 
constant. Note that varying the rapidities at fixed y\ — y 2 would keep the product xuxa constant, and would force a 
competition between the evolution of q(xh, A 2 ) with increasing Xh and the CGC evolution with decreasing xa- 




FIG. 3: The A<j> spectrum ([59} in two situations for the RHIC energy y/s = 200 GeV/nucleon. Fig.3a: p Tl = 3.5 GeV, 
Pt 2 = 2 GeV, j/i =3.5 and 2/2 is varied from 1.5 to 2.5. Fig. 3b: pr x = 5 GeV, j/i = 3.5, t)2 = 2 and pr 2 is varied 1.5 GeV to 3 GeV. 
In both cases, the correlation in azimuthal angle is suppressed as the value of xa probed in the process decreases. Varying pr 2 
at fixed 2/2 is much more efficient as the ratio Pt 2 /Qs varies over a larger range. 



In Fig. 3a, we have studied the A(f> spectrum (|59|) in the situation in which pt x = 3.5 GeV, pr 2 = 2 GeV, y\ =3.5 and 
j/2 is varied from 1.5 to 2.5. As 2/2 increases, the value of xa decreases and the suppression of the azimuthal correlation 
is more important. However the effect is quite small, because the increase of the saturation scale with decreasing x a 
is rather slow. In Fig. 3b, we investigate the situation for which px 1 — 5 GeV, 2/1 = 3.5, j/2 = 2 and pt 2 is varied 1.5 GeV 
to 3 GeV. As pr 2 decreases, it gets closer to the saturation scale Q s (which also slightly increases as xa decreases), 
and the suppression of the azimuthal correlation increases. Varying px 2 at fixed 2/2 allows to probe the ratio pt 2 /Qs 
over a larger range, so the effect is much bigger than when varying 2/2 at fixed pt 2 ■ 

Experimental measurements of two-particle correlations in azimuthal angle have been performed at RHIC in d—Au 
collisions [23] by the PHENIX and STAR collaborations. Our predictions for the fully differential cross section are not 
directly comparable with the data. One would have to carry out a number of integrations over the kinematic variables, 
while properly taking into account the kinematic cuts applied by the experiments for the different measurements, but 
this goes beyond the scope of this work. Nevertheless the exploratory measurements of STAR with tt° at forward 
rapidity and charged hadrons at mid rapidity are qualitatively consistent with a suppression of the back-to-back peak 
with respect to p— p collisions. By contrast, the measurements of PHENIX do not show any evidence of a suppression 
of the back- to-back peak, but they probe values of xa which are bigger than 0.01. It may very well be that the CGC 
picture breaks down for values of xa bigger than 0.01, and it justifies our choice not to start the small— evolution 
at a higher value. 

In Fig. 4, the A<3> spectrum (|3T)f is plotted for the same situations as in Fig. 3, but with the LHC heavy-ion energy 
^=5.5 TeV/nucleon. Assuming similar possibilities for the ALICE detector, compared to RHIC detectors, the final- 
state kinematics are unchanged, and as a result the values of xa probed in the process are much smaller at the LHC 
(typically xa ~ 5.10 -5 ,) compared to RHIC (typically xa^10~ 3 ). One sees that in both cases, the azimuthal angle 
decorrelation behaves as a function of 2/2 and pr 2 as in Fig. 3, but as indicated by the vertical scale, the spectrum is 
globally more suppressed and the peak is also slighlty broader. Let us warn that those conclusions are only qualitative, 
as our calculation is really only suited for RHIC where the q—>qg process is predominant. The values of Xh probed 
at the LHC (typically x^~0.02 compared to Xh^0.5 at RHIC) are such that the gluon-initiated processes g^qq and 
g— >gg (not included in our calculation) will dominate the cross-section. 



V. CONCLUSIONS 



Let us summarize our main results. We computed forward inclusive dijet production qT—fqgX in the scattering of 
a quark off a Color Glass Condensate. The two-particle spectrum ([24l was expressed in terms of correlators of Wilson 
lines. With the Gaussian CGC wavefunction (f3"2"| , we could compute the correlators in terms of a single function, 
which in practice is obtained (in the large— N c limit) by solving the BK equation f3"§|) with the MV initial condition 
(|4T)|) . We applied our expression |4T|) to the process hT — ► h\h%X, the inclusive production of two particles hi and 
/12 at forward rapidities, in the direction of the dilute hadron h. As an application of formula (f54|) . we studied the 
azimuthal angle correlation in d~Au collisions. While our results for the fully differential spectrum (I59|) are not yet 
comparable with the data, we obtain a qualitative agreement. 



pr t = 3.5 GcV, pt 2 — 2 GeV, yi = 3.5 PTi=5GeV, yi — 3.5, y% = 2 
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FIG. 4: The Ac/> spectrum (|59[l in the two situations studied in Fig. 3, but for the LHC energy -^5 = 5. 5 TeV/nucleon, resulting 
in probing much smaller values of xa- In both cases, the correlation in azimuthal angle varies as a function of yi and pr 2 as in 
Fig. 3, but globally the azimuthal correlation is more suppressed (see the vertical axis) and the peak is less pronounced. 



However, if the CGC discovery at RHIC is to be promoted to the same level than that of the quark-gluon plasma 
[23 | . then the qualitative agreements should be made quantitative. This has been done in the case of single particle 
production at forward rapidities in d— Au collisions (25j |. For azimuthal correlations, this paper represents a first step, 
but more efforts are required, on both the experimental and theoretical sides. Other d— Au runs at RHIC in the future 
would certainly be of interest, for instance the accessible range in rapidity could be improved. 

The present calculation takes into account the effects of the non-linear QCD evolution at small— x. However, since 
the validity of the CGC picture requires x < 0.01, with the RHIC energy the range probed in x is somewhat limited 
(down to 10~ 3 ), and the evolution might be only tested at the LHC. In this case the processes gT — > qqX and 
gT^ggX for dijet production should also be included in the calculation. At the level of formula JH]), the expressions 
of [26j ] could be useful, but to proceed further in the derivation, one needs to compute an 8-point function (for the gg 
final state), which we leave for future work. 

Finally, the final-state configuration studied in this paper requires both particles to be produced at forward ra- 
pidities, in order to avoid a large rapidity interval between them. The situation considered in Q with one particle 
produced at forward rapidity and the other at mid-rapidity calls for the inclusion of other small— x QCD effects in 
the BFKL framework [lj|. In the context of azimuthal angle correlations, these effects have been the focus of de- 
voted studies [13, however combining them with the CGC evolution included here is still an open and interesting 
problem. 
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